h expansion for the periodic orbit quantization by harmonic inversion 
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Semiclassical spectra beyond the Gutzwiller and Berry- 
Tabor approximation for chaotic and regular systems, respec- 
tively, are obtained by harmonic inversion of the h expansion 
of the periodic orbit signal. The method is illustrated for 
the circle billiard, where the semiclassical error is reduced by 
one to several orders of magnitude with respect to the lowest 
order approximation used previously. 

PACS numbers: 03.65.Sq, 05.45.+b 



Semiclassical spectra can be obtained for both regular 
and chaotic systems in terms of the periodic orbits of the 
system. For chaotic dynamics the semiclassical trace for- 



, and for integrable 
is well known to be 



mula was derived by Gutzwiller | 
systems the Berry- Tabor formula ^_ 
precisely equivalent to the EBK torus quantization [Q. 
However, the semiclassical trace formulae are exact only 
in exceptional cases, e.g., the geodesic motion on the con- 
stant negative curvature surface. In general they are just 
the leading order terms of an infinite series in powers 
of the Planck constant and the accuracy of semiclassi- 
cal quantization is still an object of intense investigation 
. Methods for the calculation of the higher order pe- 
riodic orbit contributions were developed in || 0. How- 
ever, the h expansion of the periodic orbit sum does not 
solve the general problem of the construction of the ana- 
lytic continuation of the trace formula. The semiclassical 
trace formula usually does not converge in the physically 
interesting region even when only the leading order terms 
in h are considered, and special techniques are necessary 
to overcome the convergence problems ]H|-p^[. Up to 
now the h expansion for periodic orbit quantization is re- 
stricted to systems with known symbolic dynamics, like 
the three disk scattering problem, where cycle expansion 
techniques can be applied |p|Jl0||. 

Recently the harmonic inversion technique Jl4|,[l5| was 
proposed as a universal method for periodic orbit quan- 
tization Jlj],[l7]], which allows the analytic continuation 
of the non-convergent periodic orbit sum to the region 
where the semiclassical eigenvalues and resonances are 
located. The power of this method was demonstrated by 
its wide applicability to open and bound systems with 
both regular and chaotic classical dynamics. However, 
the method was restricted to the conventional lowest or- 
der h approximation of the periodic orbit sum, i.e., it 
cannot be applied straightforwardly to the Ti expansion 
of the periodic orbit sum. In this Paper we overcome 



these problems and extend the method of periodic orbit 
quantization by harmonic inversion to the analysis of the 
h expansion of the periodic orbit sum. When applied 
to the circle billiard, as a first example, the accuracy of 
semiclassical eigenvalues is improved by at least one to 
several orders of magnitude. The method can be applied 
to a large variety of systems, i.e., it is not restricted to 
problems which can be solved with cycle expansion tech- 
niques. 

As previously [|l6| we consider systems with a scaling 
property, i.e., where the shape of periodic orbits does not 
depend on the scaling parameter, w, and the classical 
action S po scales as 



(1) 



The scaling parameter plays the role of an inverse effec- 
tive Planck constant, i.e., w = h~g , and the K expansion 
of the periodic orbit sum can therefore be written as a 
power series in w^ 1 . The semiclassical spectrum is given 
by 



q(w) = --lmg(w) , 
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the fluctuating part of the semiclassical response func- 
n „„ +u„ ^th 



tion. 



The ApJ are the complex amplitudes of the n 



order periodic orbit contributions including phase infor- 
mation from the Maslov indices. Usually the zeroth or- 
der contributions A p °o are considered only. The Fourier 
transform of the principal periodic orbit sum 



C (s) = ±- [ °° g (w)e- isw dw 
= ^4°>( S - Spo ) 
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is adjusted by application of the harmonic inversion tech- 
nique jl6|,[l7j to the functional form of the exact quantum 
expression 
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with {wk, dk} the eigenvalues and multiplicities. The fre- 
quencies Wk,o obtained by harmonic inversion of Eq. ^ 
are the zeroth order h approximation to the semiclassical 
eigenvalues. We will now demonstrate how the higher or- 
der correction terms to the semiclassical eigenvalues can 
be extracted from the periodic orbit sum (§). We first 
remark that the asymptotic expansion (^|) of the semi- 
classical response function suffers, for n > 1, from the 
singularities at w = 0, and it is therefore not appro- 
priate to harmonically invert the Fourier transform of 
(||), although the Fourier transform formally exists. This 
means that the method of periodic orbit quantization by 
harmonic inversion cannot straightforwardly be extended 
to the h expansion of the periodic orbit sum. Instead we 
will calculate the correction terms to the semiclassical 
eigenvalues separately, order by order, as described in 
the following. 

Let us assume that the (n— f) st order approximations 
Wk, n -i to the semiclassical eigenvalues are already ob- 
tained and the Wk, n are to be calculated. The difference 
between the two subsequent approximations to the quan- 
tum mechanical response function reads 



,(«0 



dk 



t-f \ W - W k ,n + it W - Wk.n-1 + ^ 
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with Wk, n = (wk.n + i0fc, n _i)/2 and Aw k , n = Wk, n - 
Wk, n -i- Integration of Eq. ^ and multiplication by w n 
yields 

g n (w)=w g n {w)dw = y = — , (7) 

J <-? w - Wk,n + ie 

which has the functional form of a quantum mechanical 
response function but with residues proportional to the 
n th order corrections Awk, n to the semiclassical eigenval- 
ues. The semiclassical approximation to Eq. |?] is obtained 
from the term g n (w) in the periodic orbit sum (||) by in- 
tegration and multiplication by w n , i.e. 



(8) 



Q n {w) = w n J g n (w)dw 
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We can now Fourier transform both Eqs. [?] and |§|, and 
obtain (n > 1) 



C„(s) = i / + G n (w)e-* ws dw 

= iJ2dk(w k ) n Aw kin e- iw * s 
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Eqs. I) and |l^ are the main result of this Paper. They 
imply that the h expansion of the semiclassical eigen- 
values can be obtained, order by order, by harmonic in- 
version (h.i.) of the periodic orbit signal in Eq. |l^ to 
the functional form of Eq. ^. The frequencies of the pe- 
riodic orbit signal ( |To| ) are the semiclassical eigenvalues 
Wk- Note that the accuracy of the semiclassical eigenval- 
ues does not necessarily increase with increasing order n. 
We indicate this in Eq. || by omitting the index n at the 
eigenvalues w k - The corrections Awk, n to the eigenvalues 
are obtained from the amplitudes, dk(wk) n Attjt lfl , of the 
periodic orbit signal. 

The method requires as input the periodic orbits of the 
classical system up to a maximum period (scaled action) , 



on), 
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•Smax, determined by the average density of states |16| 
The amplitudes Ap6 are obtained from Gutzwiller's trace 
formula G|,B and the Berry- Tabor formula || for chaotic 
and regular systems, respectively. For the next order cor- 
rection Apo explicit formulae were derived by Gaspard 
and Alonso for chaotic systems with smooth potentials 
H and in Refs. [p|Jlo[| for billiards. With appropriate 
modifications [|l8| the formulae can be used for regular 
systems as well. 

We now demonstrate the fi expansion of the periodic 
orbit sum for the example of the circle billiard. We choose 
this system mainly for the sake of simplicity, since all the 
periodic orbits and the relevant physical quantities can 
be obtained analytically. It will be evident that the pro- 
cedure works equally well with more complex systems 
where periodic orbits have to be searched numerically. 
Furthermore this system has served recently as a show- 
piece example for solving the fundamental problem of 
reducing the number of orbits required for periodic orbit 
quantization Q . The exact quantum mechanical eigen- 
values E = h k 2 /2M of the circle billiard are given as 
zeros of Bessel functions J\ m \(kR) — 0, where m is the 
angular momentum quantum number and R, the radius 
of the circle. In the following we choose R = 1. The 
lowest order semiclassical eigenvalues can be obtained by 
an EBK torus quantization resulting in the quantization 
condition |5| 



kRy/l — (m/kR) 2 — \m\ arccos 
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with m = 0, ±1,±2, ... being the angular momentum 
quantum number and n = 0, 1, 2, . . . the radial quantum 
number. States with angular momentum quantum num- 
ber m/0 are twofold degenerate (dk — 2). 

For billiard systems the scaling parameter is the ab- 
solute value of the wave vector, w = k = and 
the action is proportional to the length of the orbit, 
SpQ — 7ik£pQ. The periodic orbits of the circle billiard 
are those orbits for which the angle between two bounces 
is a rational multiple of 2n, i.e., the periods £ po are ob- 
tained from the condition 
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2m r sin 7 , 
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with 7 = 7rm^/TO r , = 1,2,... the number of turns of 
the orbit around the origin, and m r = 2m^,, Ira^ + 1, . . . 
the number of reflections at the boundary of the circle. 
Periodic orbits with m r ^ 2m^ can be traversed in two 
directions and thus have multiplicity 2. As mentioned 
before the calculation of the zeroth order amplitudes „4po' ■— . 
in Eq. || depends on whether the classical dynamics is 
regular or chaotic. For the circle billiard with regular 
dynamics we start from the Berry- Tabor formula || and 
obtain 
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where /i po = 3m. r is the Maslov index. [Note that the 
factor w~ n in Eq. || must be replaced by w^™" 1 / 2 ) for 
the regular circle billiard.] For the calculation of the first 
order periodic orbit contribution gi (w) in Eq. [3] we adopt 
the method of Alonso and Gaspard |9| and obtain the first 
order periodic orbit amplitudes 



5-2sin 7 e - i( fp po -f) 



sin 7 



(14) 



A detailed derivation of Eq. [14| will be given elsewhere 

f. With the periodic orbit amplitudes, Eqs. [l3] and 
at hand we have all the ingredients necessary for the 
harmonic inversion of the zeroth and first order periodic 
orbit signal. For the technical details of the harmonic in- 
version technique see Refs. Jl4|,|l5],[l7|] . We considered pe- 
riodic orbits up to maximum length £ max = 200, which 
was sufficient to resolve the low lying states, despite a 
few near degeneracies. The zeroth order semiclassical 
approximations k 1 -- ^ to the eigenvalues are obtained by 
harmonic inversion of the signal Cq(s) (Eq. and are 
presented in Table [|. They agree (despite the near degen- 
erate states at k 11.05 marked by asterisks) within the 
numerical accuracy with the results of the torus quanti- 
zation, Eq. [ll] (see eigenvalues /c EBK in Table ||). How- 
ever, the semiclassical eigenvalues deviate significantly, 
especially for states with low radial quantum numbers n, 
from the exact quantum mechanical eigenvalues A: cx in 
Table |. 

The first order corrections to the semiclassical eigenval- 
ues are obtained by harmonic inversion of the peri- 
odic orbit signal Ci(s) (Eq. |To|). The resulting spectrum, 
i.e., the integrated differences of the density of states 
J Ag(k)dk are shown in Fig. |l|. The squares mark the 
spectrum for Ag(k) — g^(k) — g(°\k) obtained from 
the harmonic inversion of the signal Ci(s). For compar- 
ison the crosses present the same spectrum but for the 
difference Ag(k) — g cx (k) — g EBK (k) between the exact 
quantum mechanical and the EBK-spectrum. The devia- 
tions between the peak heights exhibit the contributions 
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FIG. 1. Integrated difference of the density of states, 
J Ag(k)dk, for the circle billiard with radius R = 1. Crosses: 
Ag(k) = g cx {k)-g EBK (k). Squares: Ag(k) = g m (k)-g m {k) 
obtained from the h expansion of the periodic orbit signal. 



of terms of the h expansion series beyond the first order 
approximation. 

The peak heights of the levels in Fig. |l| (solid lines and 
crosses) are, up to a multiplicity factor for the degenerate 
states, the shifts Ak between the zeroth and first order 
semiclassical approximations to the eigenvalues k. The 

TABLE I. The 20 lowest eigenstates of the circle billiard 
with radius R = 1. n,m: Radial and angular momentum 
quantum numbers; fc EBK : Results from EBK-quantization; 
k 1 " ': Eigenvalues obtained by harmonic inversion of the pe- 
riodic orbit signal without Ti corrections (nearly degenerate 
states marked by asterisks are not fully resolved); k^: Eigen- 
values obtained by harmonic inversion of the periodic orbit 
signal including h correction; fc cx : Exact eigenvalues, i.e., ze- 
ros of the Bessel functions J m (kR) = 0. 



n 


m 


fc EBK 














2.356194 


2.356187 


2.409239 


2.404826 





1 


3.794440 


3.794430 


3.834226 


3.831706 





2 


5.100386 


5.100379 


5.138108 


5.135622 


1 





5.497787 


5.497782 


5.520501 


5.520078 





3 


6.345186 


6.345180 


6.382687 


6.380162 


1 


1 


6.997002 


6.996999 


7.015857 


7.015587 





4 


7.553060 


7.553053 


7.590944 


7.588342 


1 


2 


8.400144 


8.400140 


8.417477 


8.417244 


2 





8.639380 


8.639370 


8.653839 


8.653728 





5 


8.735670 


8.735652 


8.774088 


8.771484 


1 


3 


9.744628 


9.744619 


9.761243 


9.761023 





6 


9.899671 


9.899663 


9.938844 


9.936110 


2 


1 


10.160928 


10.160925 


10.173526 


10.173468 


1 


4 


11.048664 


11.048966* 


11.077169* 


11.064709 





7 


11.049268 


11.048966* 


11.077169* 


11.086370 


2 


2 


11.608251 


11.608248 


11.619883 


11.619841 


3 





11.780972 


11.780968 


11.791546 


11.791534 





8 


12.187316 


12.187318 


12.228099 


12.225092 


1 


5 


12.322723 


12.322717 


12.338791 


12.338604 


2 


3 


13.004166 


13.004163 


13.015235 


13.015201 
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first order eigenvalues k^ = fcW + Ak are presented in 
Table ffl, and are in excellent agreement with the exact 
eigenvalues fc ox . An appropriate measure for the accu- 
racy of semiclassical eigenvalues is the deviation from 
the exact quantum eigenvalues in units of the average 
level spacings, (Afc) av = l/g(k). Fig. || presents the 
semiclassical error in units of the average level spacings 
(Afc) av « 4/fc for the zeroth order (diamonds) and first 
order (crosses) approximations to the eigenvalues. States 
are labeled by the radial and angular momentum quan- 
tum numbers (n, m). In the zeroth order approximation 
the semiclassical error for the low lying states is about 
3 to 10 percent of the mean level spacing. This error 
is reduced in the first order approximation by at least 
one order of magnitude for the least semiclassical states 
with radial quantum number n = 0. The accuracy of 
states with n > 1 is improved by two or more orders of 
magnitude. 
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FIG. 2. Semiclassical error |fc (0) - fc cx | (diamonds) and 
fc' 1 ' — fe cx | (crosses) in units of the average level spacing 
(Afc} av w 4/fe. States are labeled by quantum numbers 
(n, m). 

As mentioned above the general technique developed in 
this Paper is not restricted to the circle billiard but can in 
general be applied to the whole variety of systems which 
can be quantized semiclassically by harmonic inversion of 
the periodic orbit sum. While for the circle billiard the 
periodic orbit parameters can be calculated analytically 
the orbits must be obtained from a numerical periodic 
orbit search in general. However, no additional periodic 
orbits need to be searched for the fi expansion of the 
periodic orbit sum, i.e., it is sufficient to calculate the 
amplitudes in Eq. || for the given set of orbits as described 
in Refs. 

In conclusion, we have demonstrated that semiclassical 
spectra beyond the Gutzwiller and Berry- Tabor approx- 
imation can be obtained by harmonic inversion of the 
h expansion of the periodic orbit signal. For the circle 
billiard, as a first example, the semiclassical error is re- 
duced by at least one to several orders of magnitude by 
just including the lowest order periodic orbit correction 



terms. The method proposed in this Paper opens the way 
to the calculation of high precision semiclassical eigenval- 
ues directly from periodic orbit data for both regular and 
chaotic systems. It is not restricted to bound systems but 
can be applied to open systems as well. 
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